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ABSTRACT 

We have carried out an in-depth study of the peripheral region of the molecular cloud L1204/S140, 
where the far-ultraviolet radiation and the density are relatively low. Our observations test theories of 
photon-dominated regions (PDRs) in a regime that has been little explored. Knowledge of such regions 
will also help to test theories of photoionization-regulated star formation. [Cii] 158 (im and [Oi] 63 
/jm lines are detected by ISO at all 16 positions along a 1-dimensional cut in right ascension. Emission 
from H 2 rotational transitions J = 2 —> 0 and J = 3 —> 1, at 28 and 17 ;jm, was also detected at several 
positions. The [Cn], [Oi], and H 2 intensities along the cut show much less spatial variation than do the 
rotational lines of 12 CO and other CO isotopes. The average [Cii] and [Oi] intensities and their ratio 
are consistent with models of PDRs with low far-ultraviolet radiation (Go) and density. The best-fitting 
model has Go ~ 15 and density, n ~ 10 3 cm~ 3 . Standard PDR models underpredict the intensity in the 
H 2 rotational lines by up to an order of magnitude. This problem has also been seen in bright PDRs and 
attributed to factors, such as geometry and gas-grain drift, that should be much less important in the 
regime studied here. The fact that we see the same problem in our data suggests that more fundamental 
solutions, such as higher H 2 formation rates, are needed. Also, in this regime of low density and small 
line width, the [0 1 ] line is sensitive to the radiative transfer and geometry. Using the ionization structure 
of the models, a quantitative analysis of timescales for ambipolar diffusion in the peripheral regions of 
the S140 cloud is consistent with a theory of photoionization-regulated star formation. Observations of 
[C 11 ] in other galaxies differ both from those of high Go PDRs in our galaxy and from the low Go regions 
we have studied. The extragalactic results are not easily reproduced with mixtures of high and low Go 
regions. 

Subject headings: Infrared: ISM: lines and bands — Infrared: ISM: continuum ISM: clouds 
(L1204/S140) — ISM: atoms — ISM: molecules 


1. INTRODUCTION 

Photon-dominated regions (or photodissociation re¬ 
gions, PDRs) are regions of the neutral interstellar 
medium (ISM) where far-ultraviolet (FUV) (6 eV < hv < 
13.6 eV) photons control the heating and chemical pro¬ 
cesses. They are the interface between H 11 regions and cold 
molecular cores. The physical and chemical structure of 
PDRs depends critically on the FUV intensity and the gas 
density. In the peripheral regions of the molecular cloud 


L1204/S140, the FUV intensity and the gas density are 
low, allowing tests of the models in an important regime. 
We are motivated by two primary goals: understanding 
PDRs in a regime of parameter space that has not been 
extensively studied; and understanding the role of the re¬ 
gions with relatively low FUV and density in global star 
formation, both in our own Galaxy and in other galaxies. 

A great deal of observational and theoretical effort has 
been devoted to understanding PDRs. Comprehensive 
models of PDRs have been constructed by several groups 
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(e.g., Black & Dalgarno 1977, van Dishoeck & Black 1986, 

1988, Tielens & Hollenbach 1985a, Sternberg & Dalgarno 

1989, le Bourlot et al. 1993, Kaufman et al. 1999) by solv¬ 
ing the full chemistry and heating-cooling balance in a self- 
consistent way. In PDRs with n < 10 5 cm' 3 , the most im¬ 
portant heating process for gas is photoelectric heating— 
electrons ejected from dust particles by FUV photons heat 
the gas through collisions. For dust, direct absorption of 
FUV photons is the primary heating mechanism. Far- 
infrared continuum emission is the major cooling process 
for dust, and the far-infrared lines of [Cii] 2 P 3/2 —> 2 Pi /2 
at 157.7409 pm (hereafter [C n]) and [Oi] 3 Pi —» 3 P 2 at 
63.183705 pm (hereafter [O i]) are the most important for 
gas. Therefore, the [C n] and [O i] lines, along with H 2 ro¬ 
tational emission, are the most important tracers of PDRs. 

Most previous work has focused on bright PDRs very 
close to hot OB stars, e.g., the Orion bar (Tielens & 
Hollenbach 1985b; Jansen et al. 1995; Hogerheijde et al. 
1995; Tauber et al. 1994), the NGC 2023 PDR (Steiman- 
Cameron et al. 1997, Draine & Bertoldi 1996), and the 
S140 PDR (Emery et al. 1996; Timmermann et al. 1996; 
Spaans & van Dishoeck 1997). Other recent studies in¬ 
clude those by Liseau et al. (1999), and the field has been 
reviewed by Draine & Bertoldi (1999) and by Hollenbach 
& Tielens (1999). These regions have a far-ultraviolet in¬ 
tensity of Go > 10 3 and a density higher than 10 4 cm -3 , 
where Go is the enhancement factor relative to the stan¬ 
dard interstellar radiation field as given by Habing (1968). 

There has been very little exploration of the physics 
of PDRs with modest far-ultraviolet fields and densities, 
conditions likely to prevail over most of the surface of 
molecular clouds in our Galaxy. Federman et al. (1995) 
and van Dishoeck & Black (1988) have studied PDRs in 
diffuse clouds ( Ay ~ 1 mag) and translucent clouds ( 
n < 1000 cm -3 , Go < 17 and Ay < 5 mag). Regions 
with high densities and moderate UV fields ( n < 5 x 10 4 
cm -3 , Go = 900) have also been studied in some de¬ 
tail (Jansen et al. 1995), and Kemper et al. (1999) have 
used submillimeter and far-infrared observations to probe 
a reflection nebula with n ~ 5000 cm -3 and Go ~ 200. 
In this paper, we explore the critical intermediate regime 
where n = n(H) + 2n(H2) ~ 500 — 5000 cm' 3 and 
Go ~ 10 — 60. The Infrared Space Observatory (ISO 2 ) 
provided a unique opportunity to observe low-brightness 
extended [C 11 ], [0 1 ], and H 2 . We used this capability to 
study the intermediate regime. 

It is also important to understand the role of regions 
with modest extinction and column density in star forma¬ 
tion. Regions like those we are studying include most of 
the mass in the interstellar medium (Hollenbach & Tielens 
1995), but their role in star formation is poorly known. 
Based on the Jeans criterion, most molecular clouds in 
the Galaxy are not sufficiently supported by thermal en¬ 
ergy and therefore should collapse under gravity to form 
stars. Such widespread collapse, however, would lead to 
a Galactic star formation rate hundreds of times higher 
than what has been observed (Zuckerman & Palmer 1974, 
Evans 1991). The observed low star formation rate seems 
to indicate that most parts of most molecular clouds are 
“sterile” (Evans 1999). Magnetic fields and turbulence are 


generally considered to play an important role in support¬ 
ing molecular clouds and preventing or slowing collapse. 
However, there is no widely accepted theory on the mech¬ 
anism of magnetic and turbulent support of molecular 
clouds. Recently, Elmegreen (2000) has argued that star 
formation does in fact happen within a few cloud cross¬ 
ing times, removing the need for cloud support. Pringle, 
Allen, & Lubow (2001) have pointed out that such a pic¬ 
ture strains methods of cloud formation, and they con¬ 
clude that visible clouds would have to form out of “dark” 
molecular matter. These scenarios beg the question: what 
prevents star formation in the great majority (~ 98%) of 
molecular gas? 

McKee (1989; Bertoldi & McKee 1996, 1997) proposed a 
mechanism of photoionization-regulated star formation to 
explain the low star formation rate in the Galaxy by com¬ 
bining two ideas: first, that magnetic fields support molec¬ 
ular clouds; and second, that PDRs occupy a large fraction 
of molecular gas. The timescale for ambipolar diffusion is 
proportional to the ionization fraction (x(e) = n(e)/n) 
and the FUV photons and cosmic rays are the sources of 
ionization. Most molecular gas is sterile because it resides 
in a PDR, where the ionization is high enough to slow am¬ 
bipolar diffusion. Only highly shielded regions are likely 
to form stars. In addition, newly formed stars inject en¬ 
ergy into the clouds, replenishing turbulence and slowing 
star formation in the rest of the cloud. In this picture, 
molecular clouds reach dynamic equilibrium when Ay ~ 8 
mag. By focusing on a peripheral region, we test the con¬ 
ditions in these regions, which should not be forming stars 
according to the theory. 

We have chosen a 1-D positional cut—we call it the S140 
cut—in the peripheral region of the molecular cloud L1204 
(Figure 1) for our PDR and chemistry study. S140, at dis¬ 
tance of ~ 910 pc (Crampton & Fisher 1974), is an Hn 
region associated with the molecular cloud L1204. In the 
rest of the paper, we will use S140 to refer to the whole 
cloud. The BOV star HD 211880 illuminates this cloud 
from the southwest side to create the visible Hu region, 
an ionization front, and a bright PDR that has been stud¬ 
ied extensively. For a picture that shows the CO emission 
in relation to the ionization front and exciting star, see 
Fig. 1 of Zhou et al. (1994). Northeast of the bright PDR 
lies the most prominent dense core (hereafter referred to 
as the S140 dense core) associated with the cloud. Nu¬ 
merous studies have been carried out on this cloud, e.g., 
Blair et al. (1978), Tafalla, Bachiller, & Martin-Pintado 
(1993), Plume, Jaffe, & Keene (1994), Emery et al. (1996), 
and Timmermann et al. (1996). The FUV radiation at 
the ionization front provided by HD 211880 is about 150 
times Habing’s mean interstellar radiation field (Keene et 
al. 1985; Habing 1968). The S140 cut, specified in Table 1, 
is far enough (~ 20') from HD211880 to have a low FUV 
radiation field (Go ~ 15 — 60, §4.1, 4.3) and modest visual 
extinction (Ay ~ 15, §4.2). 


2. OBSERVATIONS & DATA REDUCTION 


2 ISO is an ESA project with instruments funded by ESA Member States (especially the PI countries: France, Germany, the Netherlands and 
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We have observed the S140 cut with the Caltech Sub¬ 
millimeter Observatory (CSO 3 ) and ISO. From the CSO, 
we collected data on transitions of 12 CO and other CO 
isotopes, the neutral carbon [C I] 3 P 2 —> 3 Pi transition 
(hereafter [Ci]), and the HCO + J = 3—>2 line (see details 
in § 2.1). We use these lines to determine column density 
of the inner part of the cloud and to constrain the density. 
We used ISO to observe the [C n] and [O i] lines along the 
cut (see details in § 2.2), which will help us to determine 
the Go and density of the cloud, as well as the incident 
FUV radiation intensity on the clouds. Finally, we ob¬ 
served rotational transitions of H 2 with ISO at selected 
positions (§ 2.2). 

Different off-cloud positions have been used for the ISO 
and CSO observations, based on the IRAS 100 /xm map 
and the CO emission (Table 1). The ISO Off position 
had to be chosen farther away from the S140 cut than 
the CSO Off position, because of the greater extent of the 
far-infrared continuum. The CSO Off position was used in 
position-switched observations, but the ISO Off position 
was observed in the same way as the on-cloud positions 
but with longer integrations, in order to measure the [C 11 ] 
and [0 1 ] lines from the general background and establish 
appropriate baseline intensities. 

2.1. CSO Observations 

Various molecular lines and the [C 1 ] line were observed 
using the CSO, with parameters listed in Table 2. Most 
of the data were taken with a re-imaging device installed 
on the telescope so that the effective aperture was 1-2 m, 
depending on the frequency and the run (Plume & Jaffe 
1995). The beam size with the re-imager installed, listed 
in Table 2 for each run, was 2 , -3' and is referred to as 
the “big beam”. Beam sizes of the big beam in Table 2 
are accurate to better than 10%. Since our object is ex¬ 
tended, the exact beam size is not critical. The [C 1 ] line 
was observed with the full dish (when the re-imager was 
unavailable), in the mode of “On-The-Fly Mapping”; then 
the data were smoothed to 120" to match the big beam. 
The pointing and the main-beam efficiency were obtained 
by observing the available planets and the Moon. For the 
full dish observations, the pointing uncertainty was less 
than 5"; with the re-imager, it was no more than 30". 

The CSO facility SIS receivers were used for all obser¬ 
vations. The backend was a 1024-channel acoustic-optical 
spectrometer, with a bandwidth of 50 MHz. The typi¬ 
cal single-sideband system temperature was 300 K for the 
230 GHz band, 1300 K for the 345 GHz band, and 4000 
K for the 490 GHz band. Observations were done with 
position switching between the on- and the off-source po¬ 
sitions. The data are presented on the T^ scale. The rms 
noise is not uniform for the CSO data, because of differ¬ 
ent integration times, weather conditions, and instrument 
performance. At some positions, for some less abundant 
species, data were not taken because of extremely weak 
emission predicted from the more abundant isotopes and 
the usual abundance ratios. Repeated observations showed 
that uncertainties of the results are less than 30%. 

2.2. ISO Observations 


Parameters for the ISO observations are in Table 3. In¬ 
formation about ISO can be found in Kessler et al. (1996). 
We used the Long-Wavelength Spectrometer (LWS; Clegg 
et al. 1996; Swinyard et al. 1996) to observe the [C 11 ] and 
[Oi] lines in the LWS02 mode at 16 positions along the 
S140 cut. For each line at every position, multiple scans 
were obtained; each scan had 19 samples covering a width 
of 5 spectral resolution elements. Each sample had an in¬ 
tegration time of 0.4 sec. Particle hits caused a number 
of bad scans, resulting in a significant reduction of usable 
data. Three selected positions from the S140 cut were ob¬ 
served in the [C 11 ] and [0 1 ] lines with longer integrations. 
The ISO Off position of S140 was also observed in the [C 11 ] 
and the [0 1 ] lines. The data presented here have been pro¬ 
cessed through the Off-Line Processor Version 8.7 (OLP 
8.7), and the line fluxes were extracted by using the ISO 
Spectral Analysis Package Version 1.6a (ISAP 1.6a), cal¬ 
culated under the assumption that the sources were point 
sources. Processing consisted of removing bad scans and 
obvious glitches and, in a few cases, DC offsetting the con¬ 
tinuum to be closer to that of the other scans. The points 
were then averaged together using 2.5cr clipping and fit¬ 
ted with a line and baseline fit. The quoted errors are 
those given by the fit. Because the standard calibration 
assumes a point source and, in the region we mapped, the 
line emission is more extended than the beam, we multi¬ 
plied these fits by 0.84 for the 63 um [Oi] line and 0.59 
for the 157 um [Cn] line, using Table 4.10 of Gry et al. 
(2001). We then converted to surface brightness, using 
beam solid angles listed in Table 4.14 of Gry et al. (2001). 
The relevant detectors are SW3 ([0 1 ]) and LW4 ([C 11 ]), 
resulting in beam solid angles of 1.30 x 10 7 for [Oi] and 
6.65 x 10 -8 for [C 11 ]. The [Oi] solid angle is smaller than 
the nominal ISO beamsize of H6 but considerably larger 
than the diffraction limit. The [Cn] beamsize is slightly 
larger than the diffraction limit at 158 /nn. The point 
source fluxes were multiplied by the following factors to get 
surface brightness in ergs s _1 cnr 2 sr _1 : [Oi] 6.46 x 10 13 
; [C 11 ] 8.87 x 10 13 . In addition to our own cut through the 
cloud, we have reprocessed the LWS data taken toward 
the bright PDR near the far-infrared peak in S140 itself, 
originally taken for the program of Emery et al. (1996, 
TDT 09101820 in the ISO Data Archive) and for another 
program (TDT’s 82301120, 122, and 123) using the same 
corrections as for the data from our own program. 

Because the LWS has ten detectors, we also obtained 
undersampled far-infrared continuum scans covering from 
about 43 /tin to 175 /im from the detectors that were not 
on the lines. However, because of dark current uncertainty 
for the bands shorter than 100 /trn and the fringing prob¬ 
lem for the bands longer than 100 /trn (Swinyard et al. 
1996), we could not make accurate measurements of the 
far-infrared continuum from the ISO scans. Therefore in 
this work, we will use the IRAS measurements to constrain 
the far-infrared emission. 

The low-lying pure-rotational H 2 J = 2 —> 0 S(0) line 
at 28.218 /trn and the J = 3 —■» 1 S(l) line at 17.035 
/trn were observed with the ISO-SWS in the AOT02 mode 
(de Graauw et al. 1996). Typical integration times were 
~100 s per line, in which the 12 detectors were scanned 
once over the 28.05-28.40 and 16.96-17.11 /im ranges 


3 CSO is operated by the California Institute of Technology under funding from the National Science Foundation, contract 96-15025. 
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around the lines. The J = 5 —► 3 S(3) 9.66 pm and 
J = 7 —> 5 S(5) 6.91 pm lines were measured in paral¬ 
lel with the S(0) and S(l) lines, respectively. The spectral 
resolving power for point sources is 2000 (150 km s -1 ) 
at 28 pm, 2400 (125 km s -1 ) at 17 pm, 2280 (130 km 
s -1 ) at 9.7 pm and 1550 (195 km s -1 ) at 6.9 pm. The 
SWS apertures are 20" x 27" at S(0), 14" x 27" at S(l), 
and 14" x 20" at the S(3) and S(5) lines. Since the SWS 
beam is considerably smaller than the LWS beam, three 
positions separated by 30" in RA were observed for each 
of the three LWS positions and averaged to produce an 
intensity for comparison to the LWS data. 

The H 2 lines from translucent clouds are predicted to 
be weak and close to the sensitivity limit of the SWS in¬ 
strument. At this level, noise induced by charged-particle 
impacts on the detectors plays a large role. By good for¬ 
tune, these particle hits were exceptionally low during the 
orbit in which the H 2 data for S140 were taken so that the 
quality of the data (~0.2 Jy rms at 28 pm) is compara¬ 
ble to that commonly obtained in much longer integration 
times. The H 2 S(0) and S(l) lines were detected with the 
standard pipeline reduction, but special software designed 
to handle weak signals in combination with the standard 
Interactive Analysis Package was used for improvements 
in the signal-to-noise ratio. This technique also yields a 
possible detection of the S(3) line at one position. The 
details and justification of the methods used in the soft¬ 
ware are described by Valentijn & Thi (2000) and Thi et 
al. (2001). Because the mid-infrared continuum emission 
is weak, < 3 Jy, the data do not suffer from fringing effects 
caused by an inadequate responsivity function correction. 
However, the uncertainties are dominated by the removal 
of the fluctuating dark current (see Leech et al. 2001), 
leading us to adopt a calibration uncertainty of ~ 30%, 
which is propagated in the analysis. 

3. OBSERVATIONAL RESULTS 

Observational results from the CSO and ISO are pre¬ 
sented in Tables 4-7. Spectra of 12 COJ = 2 —> 1, 
13 CO J = 2 —> 1, C 18 0 J = 2 —► 1, and the [Ci] line at 
two selected positions are shown in Figure 2. Figure 3 
shows the spectra of the [C 11 ] and the [ 01 ] lines, averaged 
over the whole S140 cut, along with the H 2 S(0) and S(l) 
lines averaged over the three LWS positions where SWS 
observations were made. 

Two velocity components of Vlsr ~ —7.5 kms -1 and 
Vlsr ~ —10.5 kms -1 are clearly seen in the CSO data 
along the S140 cut. In Table 4, results for each compo¬ 
nent are given when applicable. The line widths (FWHM) 
of the C 18 0 J = 2 —> 1 transition of both components 
are around 1.3 kms -1 , implying a Doppler parameter, 
b = 0.78, where b is the 1/e half-width of the line. The 
first component is very close to the velocity observed at the 
S140 dense core. Therefore the molecular gas with this ve¬ 
locity is likely to be the cloud extension of the S140 dense 
core. The second component, strongly peaking around 
(774", 0) (Table 4) for the [Ci] line and the lines of 12 CO 
and other CO isotopes, is the same as that of core F iden¬ 
tified in CSJ = 1—>0 and NH 3 (1,1) and (2,2) by Tafalla 
et al. (1993). This dense core, apparently active in star 
formation with an associated IRAS source and outflow, is 
centered at (954", —60") (relative to the reference point 


of this work in Table 1) and extends about 2' x 2' in NH 3 
(1,1). We distinguish the emission in this component from 
the cloud by referring to it as core F; the first component 
is referred to as the cloud. We define the mean emission 
in core F as the average of the 5 positions with offsets in a 
between 576" and 972", using only the —10.5 kms -1 com¬ 
ponent. These positions include all with obvious effects 
from core F. The mean emission from the cloud is the 
mean of emission from all positions in the —7.5 kms -1 
component added to the mean emission from the —10.5 
kms -1 component outside that range of a. We add the 
two components for comparison to the ISO observations, 
which are not spectrally resolved. 

Table 5 presents results of the HCO + J = 3 —> 2 line 
at selected positions. At (774", 0), we detected a 0.2 K 
line of HCO + J = 3 —> 2 with the re-imager. Following 
this detection, the full dish was used to map an area of 
300" x 300" centered at (774", 0). The small-beam map 
revealed that the emission is clumpy. The strongest clump 
is close to the peak of NH 3 (1,1) detected by Tafalla et al. 
(1993) (see Figure 1). In Table 5 we list the small-beam 
results at the peak position. Big-beam observations of the 
HCO + J = 3 —» 2 line at other positions did not result in 
any detection up to the standard deviations given in Ta¬ 
ble 5. 

Measurements of the [C 11 ] and [0 1 ] lines are given in Ta¬ 
ble 6 , as is the ratio, [0 1 ]/[C 11 ], and the IRAS continuum 
intensities at 60 pm and 100 pm. Line measurements from 
the combined data set of all 16 positions along the S140 cut 
(Avg(16-pos)), at the S140 PDR peak position, and at the 
ISO off position are also given in the table. The ltr uncer¬ 
tainties come from the Gaussian fits performed by ISAP. 
These uncertainties are appropriate in evaluating the relia¬ 
bility of spatial variations in line strengths and line ratios. 
In comparisons of the absolute values of these quantities 
to theoretical models, the overall calibration uncertainty 
provides a more useful assessment of the significance of 
the results. The overall uncertainty results from a com¬ 
bination of model uncertainties for the primary calibrator 
(Uranus), variations in the detector dark current, and un¬ 
certainties in the value of the system solid angle at each 
wavelength (Swinyarcl et al. 1998). Considering the vari¬ 
ous correction factors and their uncertainties, we adopt a 
value of 30% for this overall uncertainty. 

Figure 4 shows the positional variations along the S140 
cut of the 13 CO J = 2 —► 1 and [Ci] lines (of two veloc¬ 
ity components separately), the [Cn] and [Oi] lines, and 
the intensity of the IRAS 100 pm continuum. The dot¬ 
ted lines for the panels of the [C 11 ] and [ 01 ] lines mark 
the line measurements of all the data combined of the 16 
positions (Avg(16-pos) in Table 6 ). Clearly, the [Cn] and 
[ 01 ] lines and the 100 pm continuum are quite flat over 
most of the S140 cut, while the 13 CO J = 2 —> 1 line dis¬ 
plays more than a factor of 10 variation in both velocity 
components, and the [C 1 ] varies by a factor of 4. The 
12 CO J = 2-»l, 13 CO J = 3->2, and C 18 OJ = 2-»l show 
strong spatial variation similar to the 13 CO J = 2—>1 line. 
In contrast, the average [C 11 ] and [0 1 ] fluxes of the 16 po¬ 
sitions (Avg(16-pos)) agree with the measurements at each 
individual position within 40%. Thus, core F is simply not 
apparent at all in the [C 11 ], [ 01 ], or far-infrared emission. 
Therefore, we will use the [Cn] and [Oi] data averaged 
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over all 16 positions as representative of the cloud. Core F 
is also not noticeable in the more limited H 2 data (Table 

7 ). 

The average [C 11 ] and [0 1 ] intensities are 8.66 x 10 5 
and 1.35 x 10” s erg s” 1 cm -2 sr” 1 respectively, and the 
average line ratio [0 1 ]/[C 11 ] is 0.16 ± 0.04. If we subtract 
the emission at the ISO-OFF position from the mean val¬ 
ues, the intensities for [Cn] and [Oi] are 6.97 x 10” 5 and 
0.89 x 10” 5 ergs s” 1 cm -2 sr -1 , with a ratio of [0 1 ]/[C 11 ] 
of 0.13. These numbers are very different from the S140 
bright PDR, where the ratio of [0 1 ]/[C 11 ] is 1.55 ± 0.10, 
based on our reduction of the archival data, using an anal¬ 
ysis like that used for our own data. For the peripheral 
regions, [C 11 ], rather than [Oi], is the primary coolant. 

The intensities of the H 2 lines are roughly constant at 
~ 10” 5 erg s” 1 cm -2 sr -1 at the three positions (Table 
7). In addition, the excitation temperature between J — 3 
and J = 2 is also constant within uncertainties at ~ 100 
K. The constancy of intensity and temperature is rather 
remarkable because one position is toward core F and the 
other two are positions near the ends of the cut, without 
emission from lines tracing molecular gas of moderate den¬ 
sity (cf. Table 4). Clearly, the H 2 emission has the same 
pattern of constancy as the [C 11 ] and [0 1 ] lines and traces 
the gas in the PDR. 

The H 2 S(3) line is only tentatively detected at one po¬ 
sition, with an excitation temperature between J = 5 and 
J = 3 of ~ 350 K. The non-detections at the other posi¬ 
tions give a comparable upper limit on the temperature. 

4. MODELING AND DISCUSSION 

In this section, we will attempt to reproduce the obser¬ 
vations by using both published models and our own PDR 
models. The most important free parameters in the mod¬ 
eling are the FUV radiation intensity, the total extinction, 
and the volume density. Therefore, we will first constrain 
the FUV intensity (§ 4.1) and the extinction and density 
(§ 4.2). Then we will compare the predicted [C 11 ] and [0 1 ] 
line intensities with the ISO observations in § 4.3, 4.4. 

4.1. The FUV Radiation Field Estimated from the 
Far-infrared Continuum 

Because dust grains absorb most of the FUV photons in¬ 
cident on the cloud and most of that energy is re-emitted 
in the form of far-infrared continuum, the intensity of the 
FUV radiation can be inferred from the far-infrared con¬ 
tinuum. Hollenbach et al. (1991) presented model predic¬ 
tions of the far-infrared continuum intensities at 100 /im 
and 60 /im for a range of incident FUV intensity. 

From Table 6 and Figure 4, the IRAS 100 /im and 60 /im 
intensities are quite uniform along most of the S140 cut ex¬ 
cept at Positions 15-16 for 60 /im and the Positions 13-16 
for 100 /im, where the intensities are lower. Therefore we 
use the intensities averaged over positions 1-13 and the in¬ 
tensity at position 16 to estimate a range of values for the 
ultraviolet intensity. The continuum intensities at the ISO 
Off position have been subtracted to yield I (60) = 30 — 50 
MJy sr -1 and /(100) = 100 — 180 MJy sr -1 . Assuming 
a A -1 dust emissivity, the total far-infrared intensity is 
about 6 x 10” 3 ergs s” 1 cm” 2 sr -1 . 

These intensities and aA' 1 dust emissivity yield a dust 
temperature of (28 ±2) K, near the constant value of 27 K 


reached for Go < 160, when effects of transient heating of 
small grains/PAHs by single photons are included (Hollen¬ 
bach & Tielens 1995). Therefore, the 60 /im continuum 
intensity in this region is probably strongly affected by 
small grains. Based on Fig. 18 of Hollenbach et al. (1991), 
we used the 100 /an intensity to estimate the strength of 
the ultraviolet field as Go = 40 — 60. As noted by Kaufman 
et al. (1999), these values assume that the FUV impinges 
only on one surface of the cloud. If the cloud is heated 
from both sides, then the values of Go appropriate to one 
surface would be 20 - 30. The assumption here is that 
the cloud is optically thin to the far-infrared continuum. 
Given the uncertainties in the fraction of the grain heat¬ 
ing caused by FUV photons, we estimate Go in the range, 
15-60. The enhancement factor of the radiation field Ijjv 
with respect to the standard field given by Draine (1978) 
is 1.71 times smaller than Go, Iuv = Go/1.71 (Draine & 
Bertoldi 1996). 

4.2. Column Density and Density Regime 

The mean integrated intensity ofC 18 OJ = 2—>1 for 
the cloud component is 0.54 K kms , From Table 4, the 
typical excitation temperature of 12 CO is about 10 K. As¬ 
suming that the C 18 OJ = 2—>1 line is optically thin and 
in LTE, we calculate analytically the total column density 
of C 18 0 to be < 1.3 x 10 14 cm -2 . The empirical relations 
of Frerking, Langer, & Wilson (1982) on p Oph yield a 
visual extinction, Ay = 3.4 mag, and a total column den¬ 
sity of H nuclei of 5 x 10 21 cm” 2 . Assuming the cloud is 
spherical and the angular size of the cloud along the S140 
cut is roughly 15', the linear size of the depth of the cloud 
would be 1.2 x 10 19 cm (4 pc) at a distance of 910 pc. 
Therefore, the average density along the S140 cut would 
be 800 cm -3 . The assumption of LTE is not valid at such a 
low density. Models of excitation, including trapping with 
LVG codes, indicate that the C 18 0 populations at such 
low densities are far from LTE. A more self-consistent so¬ 
lution is a density around 2000 cm ” 3 and a column density 
of C 18 0 of 2.0 x 10 15 cm” 2 , yielding A v = 16. The ratio 
ofJ=2 — ltoJ = 3 — 2 lines of 13 CO is roughly consis¬ 
tent with a density of 2000 cm” 3 , and the non-detection 
of HCO + J = 3 — 2 outside core F limits the density to 
n < 10 4 cm ” 3 for typical abundances in translucent clouds 
and Tk = 10 — 20 K. Similar considerations applied to the 
position of core F (774,0) lead to estimates of Ay = 25 
mag and, using a size of 3', a mean density of 10 4 cm” 3 . 

Considering the uncertainties in these estimates, the 
likely uncertainties in column density and density are fac¬ 
tors of 2 and 3, respectively. A direct comparison between 
PDR models and our data will be given in §4.4, which 
supports these simple estimates. 

4.3. Comparison to Published PDR Models 

We first compare our results to the published grid of 
PDR models by Kaufman et al. (1999). These calcu¬ 
lations incorporated new collision rates for fine-structure 
lines and H 2 , new PAH heating and chemistry, and lower 
gas-phase abundances for oxygen and carbon (Savage & 
Sembach 1996). They include regions of low FUV and 
low density, which seem to be appropriate for our situa¬ 
tion. The calculations employed a turbulent broadening 
of 1.5 kms” 1 , about twice what we infer, but Kaufman 
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et al. argue that the results are not too sensitive to this 
parameter. 

We first consider the constraints from the two atomic 
species, [Cn] and [O i]. From Fig. 4 of Kaufman et al., the 
mean ratio of [O i]/[C n] of 0.16 ± 0.4 constrains Go < 20. 
This is at the low end of the range inferred in §4.1, even 
if we assume the cloud is heated from both sides. At this 
Go, the only solution for density is n ~ 300, lower than 
the values inferred in §4.2. For lower Go, the solutions 
for n bifurcate into lower and higher values. For example, 
at Go = 10, there are solutions for log n = 2.0 or 3.3. 
The latter are more consistent with the LVG modeling of 
C 18 0 J = 2 —> 1. The [Cn] emission by itself (Fig. 3 of 
Kaufman et al.), is less diagnostic, but roughly consistent 
with the ranges implied by the ratio. 

Are these low values of Go consistent with the 
far-infrared emission? They are substantially below the 
estimates from §4.1. Our ratio of intensities of the sum of 
the [O i] and [C n] lines to the total far-infrared intensity, 
1.7 x 10 -2 , is close to the highest values found by Kauf¬ 
man et al. (1999). In fact, the ratio of line to continuum 
of the cloud in S140 may exceed the maximum in their 
models, if indeed the cloud is heated from both sides. 

What about the [Ci] data? The mean for the cloud, 
summing the two components as for the C 18 0 data, and 
correcting for efficiency, is 1.3 x 10 -6 ergs s -1 cm -2 sr -1 . 
Comparing to Fig. 7 of Kaufman et al., this value favors 
very low Go- To get Go up to even 15 requires n ~ 100. 
For the conditions that match the other lines, the [C i] 
prediction of Kaufman et al. is too high by a factor of 
3-4. 

Given reasonable uncertainties, the mean cloud data re¬ 
quire quite low values of Go and n, with most likely values 
of Go ~ 15 and n ~ 1000. These are both lower than our 
initial guesses, even if the cloud is heated from both sides. 
For this range of conditions, the “surface temperature” of 
the PDR, plotted in Fig. 1 of Kaufman et al., will lie be¬ 
tween 100 and 200 K. This is the maximum temperature 
that should apply to the region of H 2 emission, if collisions 
dominate the H 2 excitation. The excitation temperatures 
of H 2 (Table 7) are consistent with this range of tempera¬ 
tures, but the absolute intensities are difficult to reproduce 
(§ 6 )- 

Comparison of the Kaufman et al. results with those of 
other models in the literature (e.g., Le Bourlot et al. 1993, 
Roueff, private communication) and our own models (e.g., 
Jansen et al. 1995, Spaans & van Dishoeck 1997, Spaans, 
private communication), shows good agreement for the 
[C n] intensities but variations of factors of two in the [O i] 
intensities depending on the temperature structure, radia¬ 
tive transfer treatment and geometry of the source. Thus, 
the range in inferred Go and n could be somewhat larger, 
but the conclusion that both are low holds firm. 

The region of the core F emission is clearly denser than 
that of the cloud, yet core F is not apparent in either 
far-infrared lines or continuum. The emission from the 
[C n] line may be quite independent of density up to about 
n ~ 10 4 cm -3 , in the allowed range of Go; however, the 
[O i] line should increase with density. From Fig. 4 of 
Kaufman et al. (1999), the ratio of [Oi]/[Cii] should in¬ 
crease to about 1 for n = 10 4 and G 0 in the range that 
fits our data. To keep the ratio below 0.3 would require 


Go ~ 1. That is, core F would have to be shielded from 
the ambient FUV field. Interestingly, the [C i] emission 
does see core F in the —10.5 kins -1 line. Either that com¬ 
ponent is not completely shielded or the enhanced [C i] is 
caused by some internal source that does not affect the 
far-infrared line and continuum emission. 

4.4. Comparsion to PDR-Monte Carlo Models 

To check the effects of radiative transport, different 
Doppler broadening, and metal abundance at these low 
values of Go and n, we ran our own models for the part of 
parameter space indicated in the last section. We used our 
own PDR code, an updated version of the code described 
by Jansen et al. (1995) to calculate abundances of relevant 
species as a function of depth into the cloud. We then fed 
the results into a Monte Carlo code (Hogerheijde & van 
der Tak 2000) to calculate excitation and radiative trans¬ 
port. Einstein A values of [Cn] and [Oi] were taken from 
Tielens & Hollcnbach (1985a) and Galavis et al. (1997), 
respectively. In going from the PDR code to a spherical 
Monte Carlo code, one doubles the column density, so the 
PDR models were run for Ay = 9, yielding a total ex¬ 
tinction through the cloud of 18 mag. We ran models for 
Doppler parameters from 0.5 to 2.0. The integrated inten¬ 
sities were generally insensitive to this parameter. Two 
sets of values for metal abundances were used: normal 
abundances are those listed by Jansen et al. (1995) in 
their Table 2 for the ( Oph diffuse cloud; “low metals” 
corresponds to decreasing all abundances other than C, 
N, and O by a factor of 10. The most important species 
is S, which is thought to be quite depleted in molecular 
clouds compared to translucent clouds. 

The results showed that the radiative transport of [O i] 
is indeed quite sensitive to the input parameters because 
it is very subthermally excited at these low densities, but 
quite opaque (r ~ 10 for models that fit the data well). At 
very low values of Go and n, the [O i] line predicted by the 
Monte Carlo cloud can be considerably stronger than that 
predicted by the PDR code alone, after accounting for the 
doubling of the column density. The differences seem to 
be partially in the handling of the radiative transport and 
partly in the geometry. Calculations that calculate the 
temperature self consistently in spherical geometry pro¬ 
duce even larger effects (M. Spaans, personal communica¬ 
tion) . 

Among a grid of models, the best fit was obtained for 
a model with Go = 17 and n = 1000 and low metals. 
This model reproduced the C 18 0 line strength, indicating 
that the extinction estimate from §4.2 was about right. 
The ratio of 13 CO J = 3 —> 2/ 13 CO J = 2 —» 1 was about 
50% lower than the observations in this model, while mod¬ 
els with n = 3000 gave ratios higher than observed and 
13 CO J = 2—>1 lines about a factor of 2 stronger than ob¬ 
served. The HCO + J = 3—>2 line produced by this model 
is about 1.5 times our RMS noise, still consistent with the 
observations, while the model with n = 3000 would pro¬ 
duce a 10 a detection. The HCO + abundance and thus 
the HCO + J = 3^2 line are considerably enhanced in the 
case of low metals. Taken together, these results suggest 
that n = 1000 cm -3 is about right, but perhaps a bit low. 

The best model gave a [Cn] line in good agreement 
with observations and an [O i] line about 40% too strong. 
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Higher values of Go or n greatly overproduced the [Oi] 
line. These are in rough agreement with Kaufman et ah, 
but the [C i] line predicted by our models is about a factor 
of 4 weaker than the line predicted by Kaufman et al., and 
our prediction agrees with the observations. Kaufman et 
al. note that their [C i] lines are stronger than in previous 
calculations because of inclusion of PAH chemistry in their 
models. While charge transfer from PAHs to C + to create 
C° is included in our models, the details may be different. 

We use the results of the model with Go = 17 and 
n = 1000 in the next section. 

5. TESTING PHOTOIONIZATION-REGULATED STAR 
FORMATION 

Figure 5 shows the distribution of electron fraction 
(x(e) = n(e)/n, with n = 2 n(H 2 ) + n(H)) with extinc¬ 
tion from the best fitting models from §4.4. The model 
with low metals fits the C 18 0 data better and has slightly 
lower x(e). 

The timescale for ambipolar diffusion is: 

f ad — 1-6 x 10 14 a;(e) yrs, 

where x(e) is the electron abundance (McKee et al. 1993). 
Since the lifetime of molecular clouds is roughly 30 Myrs 
(Bash, Green, & Peters 1977), if f ad exceeds this time, 
the region is effectively sterile as long as it remains mag¬ 
netically sub-critical: it will be dissipated before it can 
form stars. By setting f ad to 30 Myr, we find a threshold 
ionization fraction for star formation of x(e) ~ 2 x 10 ”'. 
In Figure 5, the threshold x(e) ~ 2 x 10”' is reached at 
Ay = 6 mag for the low metal model and is not reached 
by the end of the high metal model at Ay = 9 mag. A 
cloud with a column density less than 6-10 mag will have 
little well-shielded (x(e) < 2 x 10 ” 7 ) gas to form stars. 
Thus the ionization fraction in the region we are studying 
suggests that little star formation would be expected in 
the extended cloud, if photoionization-regulated star for¬ 
mation applies (McKee 1989). 

Based on the above analysis, the localized core F along 
the S140 cut, with Ay ~ 25 mag, has substantial gas 
that is sufficiently well shielded (x(e) < 2 x 10 ” 7 ) to allow 
star formation, which has been confirmed by the outflow 
around this core (Tafalla et al. 1993). This is consistent 
with McKee’s theory of photoionization-regulated star for¬ 
mation. In comparison, most of the molecular gas in the 
extended cloud, with Ay ~ 16 mag, has a x(e) > 2 x 10 ” 7 
and is unlikely to form stars, if subjected to FUV radiation 
from both sides of the cloud. Infrared studies to look for 
evidence of star formation in the peripheral region of S140 
(away from core F) would be able to test this prediction 
and McKee’s theory. 

Since the average extinction of molecular clouds in the 
Galaxy is about 7.5 mag (Larson 1981), according to the 
photoionization-regulated star formation model, most of 
the molecular gas will not form stars, leading to the low ob¬ 
served Galactic star formation rate. The sterility of much 
of the molecular gas is supported by the apparent absence 
of distributed populations of embedded stars in NGC 2023 
(Li, Evans, & Lada 1997). The extinction in the peripheral 
region of that cloud, estimated from 13 CO J=l —>0 emis¬ 
sion is less than 10 mag, while that in the core is about 14 
mag, just enough for shielding to decrease the ambipolar 
diffusion time. More recently, Carpenter (2000) has used 


2MASS data to set upper limits on the distributed pop¬ 
ulation in several clouds, again consistent with the result 
that star formation is largely confined to clusters forming 
in gas of higher than average extinction and density. 

Of course, a long t ad is only a necessary condition for 
the picture of photoionization-regulated star formation, 
not a sufficient one. The clouds must be magnetically 
subcritical as well. Crutcher (1999) has summarized the 
data on this topic, concluding that clouds were slightly 
supercritical, but he noted that the uncertainties allowed 
the clouds to be subcritical. Shu et al. (1999) have noted 
some corrections that make the clouds even closer to the 
critical boundary, and they suggested that observed clouds 
may be self-selected to be close to marginally critical. 

6 . THE H 2 EMISSION 

The H 2 pure-rotational lines observed with the ISO- 
SWS provide a new test of the PDR models. Early ISO- 
SWS results have shown that H 2 lines up to J = 7 —+ 5 
S(5) are readily detected in traditional high density (> 10 4 
cm” 3 ), high G 0 (> 10 3 ) PDRs such as S140 at the posi¬ 
tion of the interface with the H II region (Timmermann 
et al. 1996, see review by Wright 2000). Those data give 
excitation temperatures for J = 3 — 7 of 400-700 K, much 
higher than expected from pre-ISO PDR models (Bertoldi 
1997, Draine & Bertoldi 1999). Because of the high den¬ 
sities, the H 2 level populations are thermalized so that 
the H 2 intensities directly trace the temperature struc¬ 
ture in the warm PDR layers. To remedy these discrepan¬ 
cies, Bertoldi (1997) suggested that either the photoelec¬ 
tric heating efficiency needs to be increased or that the H 2 
formation rate on grains needs to be larger at high tem¬ 
peratures, shifting the H—>H 2 transition zone closer to the 
warm edge. Weingartner & Draine (1999) have calculated 
increased photoelectric heating efficiencies based on an en¬ 
hanced dust-to-gas ratio in the PDR due to gas-grain drift, 
and Draine & Bertoldi (2000) have shown that such mod¬ 
els can reproduce the H 2 observations of NGC 2023. Note 
that these comparisons include a factor of 5-10 enhance¬ 
ment of the models due to limb-brightening/geometry. 

In the low-density, low Go regions of S140 studied here, 
the H 2 populations should be controlled by ultraviolet 
pumping rather than collisions, at least in the outer layers 
of the cloud. Moreover, enhanced efficiencies due to gas- 
grain drift should not play a role, nor should there be any 
geometrical effects in this face-on cloud. Thus, the ISO- 
SWS data should provide a test of the H 2 excitation in a 
different PDR regime (Thi et al. 1999). For the PDR pa¬ 
rameters derived in the previous sections, n ss 1000 cm ” 3 
and Go ~ 15, standard low density PDR models such as 
those by Hollenbach et al. (1991) or le Bourlot et al. (1993) 
and our own models give H 2 excitation temperatures de¬ 
rived from the S(0) and S(l) lines of ~100-200 K, con¬ 
sistent with our ISO data. However, these models under¬ 
produce the absolute H 2 intensities by at least an order of 
magnitude. The same conclusion was reached by Habart 
et al. (2000) and Kemper et al. (1999) for somewhat denser 
PDRs exposed to low Go- Habart et al. proposed to in¬ 
crease the H 2 formation rate on grains by more than a 
factor of 5 compared with the standard H 2 formation rate 
of ~ 3 x 10 ” 17 cm 3 s ” 1 at 100 K (Jura 1975) to reproduce 
their observations. 
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If the H 2 formation rate in our S140 models is increased 
by a similar factor of 5, the comparison with observations 
is improved considerably, although the models still fall 
short by a factor of a few. Contributions from the far-side 
of the cloud could help to further reduce the discrepan¬ 
cies. Such an arbitrary, large increase of the H 2 formation 
rate is inconsistent, however, with the mechanism of the 
formation of H 2 by diffusion of physisorbed H atoms on 
silicate grains, which is effective only at very low temper¬ 
atures (Pirronello et al. 1997). A mechanism in which H 
atoms are chemically attached to the grains and in which 
the second H atom lands directly on the attached H atom, 
the so-called Eley-Rideal mechanism (Herbst 2000) is one 
possible mechanism to explain a large formation rate at 
high temperatures. Porous (carbonaceous) grains and/or 
large PAH molecules appear to be the best candidate ma¬ 
terials, but quantitative experiments and calculations on 
interstellar grain analogs at high temperatures are needed. 
In any case, an increase in the surface to volume ratio over 
the standard grain models may be needed to allow the high 
formation rates that seem to be necessary to explain our 
observations. 

An alternative mechanism to increase the H 2 (J) pop¬ 
ulations has been proposed by Spaans (1995), Joulain et 
al. (1998) and Falgarone et al. (2000)— using the energy 
released by the dissipation of turbulence or by weak (few 
km s -1 ) shocks. The observed H 2 intensities toward S140 
are within a factor of a few of those predicted by Joulain 
et al. for a cloud with n = 30 cm -3 . At densities typical 
for the periphery of S140, however, the turbulent model 
intensities should be reduced significantly. The models of 
Joulain et al. do not take ultraviolet pumping into account, 
which will enhance the populations of the higher J levels. 
Further development of these models and constraints from 
other observed species will be needed to test them fully. 

It should be noted that the strengths of the S(0) and 
S(l) lines observed in the outskirts of S140 are compara¬ 
ble to those inferred from ultraviolet observations of H 2 in 
diffuse and translucent clouds. Consider as an example the 
recent FUSE results for the translucent cloud toward HD 
73882 {Ay= 2.4 mag) by Snow et al. (2000). The observed 
column densities in J = 2, 3, and 5 translate into intensi¬ 
ties of 2.1 x 10 -6 , 1.4 x 10 -5 , and 1.7 x 10 -5 erg s _1 cm -2 
sr _1 for the S(0), S(l), and S(3) lines, respectively. The 
S(l) intensity is comparable to our observed values, but 
that of S(0) is lower and of S(3) higher. Indeed, the cor¬ 
responding excitation temperatures for the J = 2-7 levels 
of ~300 K are larger than those of ~ 100 K found here for 
the J = 2-3 levels, although they are consistent with that 
found for the J = 3-5 levels at one position. More sensitive 
searches for the higher rotational lines of H 2 are needed to 
constrain its excitation and determine the relative roles 
of ultraviolet pumping, collisions, and the H 2 formation 
mechanism in establishing the level populations. FUSE 
studies of the H/H 2 transition in well-characterized, low 
H 2 column density regions are also interesting to provide 
further constraints on the empirical H 2 formation rate in 
different environments using the analysis of Jura (1975). 

In summary, the H 2 excitation temperatures and the 
widespread nature of the emission suggest a PDR origin, 
in which the H 2 lines arise from the outer warm layers. 
However, the absolute model intensities fall short by an 


order of magnitude, unless the H 2 formation rate is en¬ 
hanced significantly. Other explanations invoked for high- 
density PDRs, such as geometry or gas-grain drift, do not 
apply to these regions. Together with the results from 
ultraviolet observations, our results suggest that our un¬ 
derstanding of the basic chemical processes involving H 2 
is still incomplete. 

7. LESSONS FOR OTHER GALAXIES 

The [C 11 ] 158 /.irn line and [0 1 ] 63 /um lines are the 
dominant coolants of the interstellar gas in galaxies. A 
significant fraction of the [Cn] emission from late-type 
galaxies arises in dense PDRs with some additional con¬ 
tributions from atomic gas not associated with molecu¬ 
lar clouds. In the Milky Way, the [C 11 ] emission fraction 
is about equally divided between PDRs associated with 
molecular clouds and neutral, purely atomic clouds (Petu- 
chowski & Bennett 1993, Mochizuki & Nakagawa 2000, 
but see Heiles 1994 for a different view). PDRs contribute 
a larger fraction of the [C 11 ] emission from galaxies with 
more active star formation (Crawford et al. 1985, Stacey 
et al. 1991) while, particularly in the outer parts of galax¬ 
ies with lower star formation rates, purely atomic gas may 
contribute to the [Cn] emission (Madden et al. 1993). 
The far-infrared continuum emission from galaxies comes 
primarily from the same PDRs, although, in the lowest 
luminosity spirals, optical emission from late-type stars 
may also contribute to the heating of dust in the molec¬ 
ular ISM (Smith et al. 1994) and something like half the 
far-infrared emission may come from purely atomic clouds 
(Mochizuki & Nakagawa 2000). It should be possible, in 
principle, to derive information about the physics of the 
neutral ISM in galaxies from the far-infrared continuum 
and fine-structure line measurements (Wolfire, Tielens & 
Hollenbach, 1990, Stacey et al. 1991, Carral et al. 1994, 
Lord et al. 1996, Malhotra et al. 2001). 

Recent studies of starburst galaxies (Luhman et al. 
1998) and of a range of lower luminosity spirals (Malho¬ 
tra et al. 1997, 2001) have made much of an apparent 
“deficit” in the [C 11 ] emission. In many of the starbursts 
and in a fraction of the spirals, the [Cii]/FIR luminos¬ 
ity ratio lies as much as a factor of 10 below the canonical 
value (3xl0 -3 ), more typical for late-type galaxies. Galac¬ 
tic photodissociation regions with high n or high Go can 
have [Cn]/FIR ratios as low or lower than those of galax¬ 
ies with the biggest [C 11 ] deficits (e.g. Genzel, Harris, & 
Stutzki 1989). The [O I] line then provides a much larger 
fraction of the cooling power (cf. the values for the S140 
cut and the S140 IRS 1 position in Table 6), and dust 
temperatures are higher. However, Malhotra et al. (2001) 
argue that the total fine structure emission ([0 1 ] plus [C 11 ]) 
is decreased substantially only if Go/n is increased. This 
combination leads to highly charged grains and inefficient 
photoelectric heating (Tielens & Hollenbach 1985a, Bakes 
& Tielens 1994). Malhotra et al. (2001) considered but 
rejected a variety of explanations, concluding that high 
values of Go/n in the emitting regions provide the best 
explanation for their sample of “normal” spirals. 

Studies of far-infrared line emission from other galax¬ 
ies have made use of the same models as adopted here, 
which assume a uniform incident radiation field and a sin¬ 
gle PDR density for the emitting material (Kaufman et al. 
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1999, Stacey et al. 1991, Wolfire et al. 1989). In addition 
to using simple models, we can gain insight into the nature 
of the regions in galaxies responsible for far-infrared fine- 
structure line emission by comparing galaxy PDR emis¬ 
sion to emission from Galactic sources with an appropriate 
range of properties. The ISO observations of fine-structure 
line emission from other galaxies take in entire giant molec¬ 
ular cloud complexes and the nuclear regions of galaxies. 
They therefore contain contributions not only from the 
high excitation regions near OB stars but also from large 
regions with lower UV fields. With our measurements of 
the cut through the outer part of the S140 cloud, a com¬ 
parison of Galactic and extragalactic fine-structure line 
emission can include a range of Galactic PDR conditions. 
Figure 6(a-c) reproduces the plots shown in Figures la, 
5, and 8 of Malhotra et al. (2001) and shows the varia¬ 
tion of Log(L([Cn])/L(F/i?)) (Fig. 6a), L([Oi])/L([Cii]) 
(Fig. 6b), and Log((L([C n]) + L([0 i]))/L(FIR)) (Fig. 
6 c) versus F 1/ (60/mr)/F ;y (100/rm) for the 45 normal, star¬ 
forming galaxies in Table 7 of Malhotra et al. (2001). Fig¬ 
ure 6 compares these data to the average values from our 
S140 cut and to typical values for high Go Galactic PDR 
sources (labeled Hu for brevity), taken from Genzel, Har¬ 
ris, & Stutzki (1989) to be L([Cn])/L(FIR) = 8.5 x 10" 5 , 
L([Oi])/L([Cii]) = 5.8, and F„(60nm)/F„(100nm) = 1.2. 
Figure 6 clearly shows that the galaxies inhabit a part of 
parameter space that is different from either the S140 cut 
or the high Go PDR sources. 

Because the galaxies generally lie between the points 
for the S140 cut and the high Go PDRs, we attempted 
to reproduce the galaxy observations by mixtures of the 
S140 cut and the high Go PDRs, indicated by the mix¬ 
ing line connecting these two points. The line represents 
parameters calculated by varying the contributions of the 
two connected Galactic components in proportion to the 
fraction of the [C n] emission contributed by each. It is 
not possible, however, to produce the observed positions 
of the galaxies in the figures with any linear combination 
of the S140 cut and the high Go PDRs. The result is 
similar if we use the emission from NGC 7023, represent¬ 
ing reflection nebulae, instead of the Hii regions. At all 
values of F„(60/rm)/i%(100/im), the mixture produces too 
much fine-structure line emission (in particular, too much 
[C n]), when compared to the galaxy sample. At a given 
value of F !y (60/rm)/F !/ (100/.tm), for example, the value of 
Log((L([Cn]) + L([Oi]))/L(F/i?)) observed for the galax¬ 
ies lies on average 0.3 below the constructed mixture. 

The results in Figure 6 show that one cannot produce 
the observed emission ratios in galaxies by simply com¬ 
bining regions like the S140 cut with regions of active 
star formation in our Galaxy. Malhotra et al. (2001) 
argued that models with high Go/n can explain individ¬ 
ual galaxies where the [CII] and [OI] is weak compared to 
the far-IR continuum. The systematically lower values of 
the galaxy ratios compared to positions along our mixing 
line, however, suggests that the conditions in heteroge¬ 
neous clouds would not lead to appropriately weak line 
fluxes at a given value of F, y (60/.im)/F 1/ (100/xm). Rather 
than focussing solely on the line to continuum ratio, it 
may be worthwhile to examine possible reasons why the 
F !y (60/im)/F !y (100/.tm) ratio might vary. If that flux ratio 
increased less rapidly along the mixing line, the data could 


be matched. Mixing very small amounts of high Go PDR 
emission has a big effect on F 1/ (60/nn)/F !y (100/.im), while 
the effect on the line emission is more modest, as indicated 
by the tick marks on the mixing line that indicate where 
90% and 50% of the [Cii] emission comes from the S140 
cut. Moderating the increase in F„(60/im)/i%(100/zm) as 
one mixes in small amounts of high Go PDR emission prop¬ 
erties could allow the mixing line to pass through most of 
the points. This would require a component of emission 
from cooler dust. Lowering the 60-100 /iin color temper¬ 
ature of the high Go PDRs from ~ 53 K to ~ 35 K (for a 
A -1 emissivity law) could account for the difference. Even 
colder dust is common in the molecular clouds near the 
center of our Galaxy (Lis et al. 2001). 

There are at least two possible causes for the cooler 
dust. The dust properties in the inner few kiloparsecs 
of galaxies, where most of the ISO emission arises, might 
differ significantly from those in nearby high Go Galac¬ 
tic PDRs. Alternatively, there might be a contribution 
from regions heated by cooler stars or from very opaque 
regions in molecular clouds. These regions would need to 
increase the fraction of the emission in the 100 /im IRAS 
band from 25% to 40% of the total far-infrared luminosity. 
Whichever reason is correct, observations of some regions 
in the inner part of our own Galaxy (GO. 68-0. 2, ISO data 
archive) find a far-infrared line to continuum ratio lower 
than for all the galaxies in the Malhotra sample, and yet 
the F, y (60/.im)/F 1/ (100/rm) is only 0.2. A mixing line con¬ 
necting GO. 68-0. 2 to the Hii region/high Go PDR’s would 
run below the bulk of the galaxies in Figures 6a and 6c. 
More complex 3-part mixtures of cold cloud edges like the 
S140 cut, cold inner galaxy clouds like GO.68 and warm 
H ii regions could presumably fill the part of parameter 
space occupied by the galaxies from the Malhotra sample. 
Understanding the properties of inner Galaxy regions like 
GO.68 may provide further insight into the conditions in 
extragalactic PDRs. 

8. SUMMARY 

We have studied the peripheral region of the molecular 
cloud L1204/S140, where the FUV radiation and density 
are relatively low, using ISO to observe [Cii], [Oi], and 
H 2 lines, and the CSO to observe the submillimeter line 
of atomic carbon and millimeter emission lines of 12 CO 
and other CO isotopes. We analyzed the results with pub¬ 
lished PDR models and our own PDR models, coupled to 
a Monte Carlo simulation code for radiative transfer. The 
important results from this work are: 

1. We detected wide-spread [Cii], [Oi], and H 2 emis¬ 
sion in the region. The [C 11 ] and [0 1 ] lines display 
much less spatial variation than [C 1 ] and the lines 
from 12 CO and other CO isotopes, which reveal the 
presence of a denser, shielded core that does not 
emit significantly in the [C 11 ] and [0 1 ] lines. 

2. The average [C 11 ] and [0 1 ] intensities of the S140 cut 
are (8.66T1.27) x lCT 5 and (1.35±0.34) x 10“ 5 ergs 
s -1 cm -2 sr -1 ; the average line ratio ([O i]/[C 11 ]) is 
0.16 ± 0.04. Therefore the [Cii] line is the domi¬ 
nant cooling line in this region. This is quite dif¬ 
ferent from what has been observed in the S140 
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bright PDR by ISO, where the ratio [O i]/[C n] is 
1.55 ±1.10. 

3. Using the model predictions in Hollenbach et al. 
(1991) and the 100 /iiri IRAS flux, the FUV intensity 
Go for the region would be 40-60, in the situation of 
single-sided FUV radiation on the cloud. If the FUV 
impinges on both sides, the inferred Go ~ 20 — 30. 
Based on C 18 0 emission and simple LVG model¬ 
ing, the density of the gas, n 1000 cm 3 , and 
Ay ~ 16. Conditions in the shielded core are differ¬ 
ent: b~1x 10 4 cm -3 , and Ay ~ 25. 

4. The PDR models (both published and our own) in¬ 
dicate that the observed [O i] intensity and the ratio 
[Cii]/[Oi] are best reproduced with Go ~ 15 and 
n ~ 1000 cm -3 . Our models, coupled to a Monte 
Carlo simulation of the radiative transfer, reproduce 
the observed C 18 0 J = 2 —*■ 1 emission for Ay ~ 18 
mag. These conditions are also consistent with the 
constraints on density from line ratios of 13 CO and 
HCO + . Our models also reproduce the observed 
[Ci] line, while the published models of Kaufman 
et al. (1999) overproduce this line. In this regime 
of n and Go, the [Oi] intensity is very difficult to 
model because it is extremely subthermal and also 
very opaque (the optical depth is 9-14). Trapping 
is very important in determining the emergent [O i] 
intensity. 

5. Quantitative analysis of the ionization structure and 
timescale for cloud collapse for the peripheral re¬ 
gions of S140 leads to conclusions consistent with a 
picture of photoionization-regulated star formation. 

6 . The detection of widespead H 2 rotational emission 
is not readily explained with standard PDR models 
and H 2 formation rates. Possible solutions include 
a substantial enhancement in the the H 2 formation 
rates or localized heating by weak shocks associated 
with turbulent decay. 

7. While the emission in [C 11 ] and far-infrared contin¬ 
uum from other galaxies lies between the charac¬ 
teristics of the emission from the peripheral region 
studied here and regions of higher Go and n, simple 
mixtures of these regions are not able to reproduce 
the characteristics of the emission from other galax¬ 
ies. 
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Fig. 1.— The S140 cut is marked as the line between two “X” symbols on the IRAS 100 pm map of the region. Darker color indicates 
stronger 100 pm emission. The open triangle is the CSO Off position, and the solid triangle marks the ISO Off. The open cross and the solid 
star are the centers of NH 3 (1,1) cores F and F 7 identified by Tafalla et al. (1993). The F 7 core has the same velocity (Vlsr ~ —7.5 kms -1 ) 
as the dense core at IRS1, while the F core is at a different velocity, Vlsr ~ —10.5 kms -1 . The star marks the position of an IRAS source 
associated with the F core. The center of the ionization front is roughly located at 22 h 17 m 30 s , 63° 01 7 ; the exciting star for the ionization 
front is at 22 h 16 m 50 s , 62° 58'. At the distance of S140, 1000" corresponds to 4.4 pc. 
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Fig. 2.— Spectra of J= 2—>1 lines of 12 CO, 13 CO, C ls O, and [Cl] 492 GHz line at two selected positions. At the position (180 // , 0), the 
line velocity (~ —7.5 km/s) is the same as that at the dense core at IRS1; thus the lines here are from the cloud extension of the dense core. 
At the position (774 // , 0) two velocity components show up, one being the same as that of the dense core, another at about —10.2 km/s. The 
second component is from a distinct clump, as evidenced by the study of Tafalla et al. 1993 
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Fig. 6 .— Ratios of line and continuum emission from galaxies in the study of Malhotra et al. (2001) compared with the same ratios in 
Galactic sources. The filled squares show the ratios for galaxies detected in both the 158 ^m [CII] line and the 63 ^m [O I] line. Where 
applicable, the open squares show the ratios for galaxies observed in [CII] but not in [O I]. The open triangles show the locations for 
galaxies detected in [CII] and observed but not detected in [O I]. In Figure 6c, the detected [CII] emission places a lower bound on the 
([Ol] + [CII])/FIR flux ratio. In this panel, the upper limit arrows end at the value corresponding to the measured [ClI ]/FIR ratio. The 
open pentagons show the limits on the ratios for galaxies observed but detected in neither [CII] or [Ol]. In each panel, the measured ratios 
for the S140 cut are shown as open circles and labeled. The measured ratios for the mean of a sample of high Go PDRs, labeled HII are 
also shown. The mixing line connecting these points is based on the fraction of [CII] emission contributed by regions like the S140 cut; the 
horizontal tickmarks along the curves indicate where such regions contribute 90% and 50% of the [CII] emission to the mixture. The vertical 
tickmarks show where regions like the S140 cut contribute 90% and 50% of the far-infrared emission. 




18 


Li et al. 


Table 1 

The ISO Positions (1950) of S140. 


Object 

S140 

Ref. Coord. 

(22 h 17 m 42 s , 63°23'15") 

End Position 1 

(1368", 0) 

End Position 2 

(-1206", 0) 

Step 

99" 

CSO Off Pos. 

(22 h 13 m 22 s , 63°23'15") 

ISO Off Pos. 

(22 h 05 m 46!7, 63°23'15") 


Table 2 


Log of CSO Observations 


Time 

Object 

Transition 


Frequency 

(GHz) 

Av a 

(kms -1 ) 

Mode b 

e b c 

(arcsec) 

TJmb^ 

June 1995 

S140 

12 CO J = 2- 

>1 

230.53797 

0.23 

RIMG 

120 

0.80 



13 CO J = 2- 

>1 

220.39870 

0.24 

RIMG 

120 

0.80 



C 18 OJ = 2- 

>1 

219.56032 

0.23 

RIMG 

120 

0.80 



13 CO J = 3- 

+ 2 

330.58800 

0.12 

RIMG 

110 

0.80 

June 1996 

S140 

[Ci] .7 = 1- 

• 0 

492.16070 

0.08 

FD 

14 

0.45 

April 1997 

S140 

HCO+J = 3- 

>2 

267.55762 

0.12 

RIMG 

160 

0.80 



HCO+J=3- 

>2 

267.55762 

0.12 

FD 

26 

0.70 


a The velocity resolution of the data. 

b The setup of the observing mode. RIMG stands for the re-imager mode and FD the full dish. 

c The beam size. 

d The main-beam efficiency. 


Table 3 

Lines Observed by ISO 


Line 


Wavelength 

(H 

E up /k a 

(K) 

E low /k h 

(K) 

^crit 

(cm -3 ) 

A 

(s- 1 ) 

Instr. 

Beam d 

[Ol] 8 Pi — > 3 

P 2 

63.183705 

227.72 

0 

4.7(5) 

8.87(-5) 

LWS 

1.3 x 10“ 7 

[Cn] 2 P 3 /2- 

H>2 Pl/2 

157.7409 

91.22 

0 

2.8(3) 

2.4(-6) 

LWS 

6.65 x 10 -8 

H 2 J = 2 -► 

0 

28.218 

509.88 

0 

5.4(1) 

2.9(-ll) 

SWS 

1.3 x 10 -8 

H 2 J = 3 -» 

1 

17.035 

1015.12 

170.48 

1.1(3) 

4.8(-10) 

SWS 

8.9 x 10 -9 

H 2 J = 5 -f 

3 

9.662 

2503.82 

1015.12 

1.9(5) 

9.8(-9) 

SWS 

6.6 x IQ" 9 


a The energy of the upper state of the transition relative the ground state in temperature. 
b The energy of the lower state of the transition relative the ground state in temperature. 
c The critical density is A/y, where 7 is the collision rate coefficient. 
d The LWS beam solid angles are taken from Gry et al. (2001). 
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Table 4 


Results of CSO Observations. 


No. 

Position 

12 C0 2—»l a 

13 CO 2—> l b 

13 CO 3—>2 b 

C 18 0 2 

->l b 

[C 

i] b 


(K) 

(K kms -1 ) 

(K km 

s -1 ) 

(K km 

s -1 ) 

(K km 

s -1 ) 

1 

(-117", 0) 


1.1 

0.6 








2 

(-18", 0) 

2.7 

1.8 

0.8 








3 

(81", 0) 

3.9 

3.4 

1.2 




0.1 




4 

(180", 0) 

6.7 

3.2 

3.9 


1.4 


0.4 


1.6 


5 

(279", 0) 

4.5 

5.4 

6.2 


1.3 

0.5 

0.4 


2.8 


6 

(378", 0) 

4.1 

5.8 

6.5 


0.8 


0.3 


3.2 

1.5 

7 

(477", 0) 

7.2 

5.1 

4.3 

2.2 

1.1 

0.9 

0.2 

0.1 

3.7 

1.5 

8 

(576", 0) 

7.2 

9.2 

6.4 

5.6 

1.5 

1.6 

0.3 

0.4 

3.7 

4.3 

9 

(675", 0) 

7.7 

9.5 

6.1 

10.6 

1.4 

2.6 

0.5 

1.4 

2.7 

3.8 

10 

(774", 0) 

8.6 

11.8 

7.3 

15.7 

2.5 

5.7 

1.2 

2.8 

2.6 

6.3 

11 

(873", 0) 

6.5 

11.3 

4.4 

14.8 

2.5 

6.2 

0.4 

2.5 


6.0 

12 

(972", 0) 

6.3 

10.2 

3.6 

8.0 

1.7 

0.8 

0.1 

1.4 


3.0 

13 

(1071", 0) 

6.6 

7.9 

3.0 

2.4 




0.2 


3.9 

14 

(1170", 0) 

5.6 

3.6 

1.3 

1.3 






1.5 

15 

(1269", 0) 

3.7 

4.1 


2.7 






1.6 

16 

(1368", 0) 

3.3 

2.4 


1.5 






1.2 

“Cloud” c 

5.6 

4.0 

4.0 

2.0 

1.6 

0.70 

0.39 

0.15 

2.9 

1.9 


a The peak of the two velocity components with Vlsr ~ —7.5 kms 1 and Vlsr ~ —10.5 kms 1 
respectively. 

b The integrated line intensities of the two velocity components in the region. 

c The mean integrated line intensities; for the —10.5 kms -1 component, the positions with 6a between 576 
and 972 were excluded. 


Table 5 

Results of the HCO + J = 3^2 Line Observations. 


Position 

FWHM Beam Area 3, 

a (rms) 


(K kms -1 ) 

(K) 


(774", 0) 

160 

0.09 0.33 

0.05 

(915", —71") b 

26 

5.75 

0.41 

(1170", 0) 

160 


0.04 

(180", 0) 

160 


0.04 


a The integrated line intensities of the two components 
are given if applicable, with the first one at Vlsr ~ —7.5 
kms -1 and the second at Vlsr ~ —10.5 kms -1 . 

b Peak position of the HCO + J = 3 —> 2 emission when 
mapped with the full dish. The HCO + J = 3 —> 2 emission 
is clumpy in the area of 300" x 300" centered at (774", 0). 
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Table 6 

Results of ISO LWS Observations. 


No. 

Position 

[Cn] 

[Oi] 

Ratio 

IRAS(60 /im) a 

IRAS(100 H a 

(arcsec) 

(10 -5 erg cm“ 

" 2 s _1 sr _1 ) 

[Oi]/[Cii] 

(MJy 

sr -1 ) 

1 

(-117,0) 

8.52±0.18 

1.05i0.06 

0.12 

57.3 

219 

2 

(-18,0) 

8.34±0.18 

0.88i0.19 

0.11 

57.2 

221 

3 

(81,0) 

8.96±0.18 

1.84A0.12 

0.21 

58.8 

220 

4 

(180,0) 

9.05±0.27 

1.26i0.19 

0.14 

66.1 

244 

5 

(279,0) 

9.05T0.27 

1.04i0.17 

0.12 

61.5 

236 

6 

(378,0) 

9.76T0.53 

1.67i0.25 

0.17 

56.2 

232 

7 

(477,0) 

7.80T0.35 

0.81i0.18 

0.10 

52.9 

226 

8 

(576,0) 

7.30T0.27 

1.09A0.21 

0.15 

52.5 

226 

9 

(675,0) 

9.58T0.27 

1.25A0.19 

0.13 

54.3 

228 

10 

(774,0) 

10.6i0.53 

1.72i0.19 

0.16 

58.8 

230 

11 

(873,0) 

8.69T0.18 

1.81A0.23 

0.21 

63.4 

227 

12 

(972,0) 

8.61T0.36 

1.68i0.32 

0.20 

63.0 

224 

13 

(1071,0) 

9.58T0.44 

1.47A0.23 

0.15 

58.3 

208 

14 

(1170,0) 

9.94T0.18 

1.55A0.11 

0.16 

51.3 

191 

15 

(1269,0) 

7.45T0.18 

1.44i0.19 

0.19 

44.8 

164 

16 

(1368,0) 

5.32T0.18 

1.03i0.21 

0.19 

39.6 

145 


Avg(16-pos) b 

8.66A1.27 

1.35i0.34 

0.16T0.04 




ISO-OFF 

1.69A0.18 

0.46A0.20 

0.27 

9.6 

47.7 


IRS1 C 

44i2 

68A3 

1.55 i 0.10 




a The IRAS fluxes at 60 /tm and 100 /jrn at the S140 cut positions and the S140 ISO off position. 

b The average line fluxes and ratios of all 16 positions. The errors are the standard deviation of the 
values at the 16 positions. 

c The average of the 4 existing measurements toward this position. The error reflects the standard 
deviation of the individual measurements. 


Table 7 

Results of ISO SWS Observations. 


No. 

Position 

H 2 J = 2 -> 0 

H 2 J = 3 —> 

1 H 2 J = 5 —> 3 

rp a 
-L 32 

T 53 a 


(arcsec) 

(io- 

" 5 erg cm -2 s" 

- 1 sr _1 ) 

(K) 

(K) 

1 

(-117,0) 

1.0A0.5 

1 .1A0.6 

< 0.6 

109 ± 30 

< 320 

“F” 

(743,0) 

0.8A0.4 

0.8T0.3 

< 0.6 

106 ± 25 

< 320 

14 

(1170,0) 

1.3A0.4 

0.9T0.4 

0.8 i 0.3 

98 i 18 

350 ± 50 


a The excitation temperature between upper J levels that yields the observed line ratio, 
assuming that the ortho-para ratio is in LTE. 



